TWO-DIMENSIONAL NONLINEAR REGIME IN THE PLASMA 
WAKEFIELD ACCELERATOR* 

A. G. Khachatryarr , S. S. Elbakian, Yerevan Physics Institute, Yerevan, Armenia 



Abstract 

The effect of nonlinearity on plasma wake wave excited 
by a relativistic cylindrical charged bunch is 
investigated. It is shown that owing to the nonlinearity 
the amplitude of wake wave gets modulated in the 
longitudinal direction. The nonlinear wavelength in the 
excited field changes in the transverse direction with the 
result that the phase front is distored and a turbulence 
developed. The nonlinear phase front distortion may be 
compensated by radial change of unperturbed plasma 
density. 

1 INTRODUCTION 

The plasma-based accelerator concepts are presently 
actively developed both theoretically and 
experimentally. This is due to the ability of plasma to 
support large acceleration rates that will reach tens of 
GeV/m, far in excess of those attained in conventional 
accelerators. Charged bunches can be accelerated and 
focused by the field of relativistic plasma waves that are 
excited by relativistic charged bunches (Plasma 
Wakefield Accelerator (PWFA)). Both the linear 
wakefield and ID nonlinear theories have been studied 
in sufficient detail (see e. g. overview in Ref. [1]). The 
allowance for finite transverse sizes of the drivers (that is 
more realistic case) and, accordingly, the transverse 
motion of plasma electrons complicate the treatment of 
the problem in the nonlinear regime. In the general case 
the analytical solution of this regime seems impossible 
and the use of numerical methods are usually required. 
Here we study the effect of nonlinearity on two- 
dimensional plasma wake waves as well as discuss the 
cause of radial steepening of the field shown in Ref. [2]. 

2 BASIC EQUATIONS 

We shall consider nonlinear steady fields excited in 
plasma by a rigid cylindrical electron bunch in the 
framework of cold hydrodynamics with immobile ions. 
Let the bunch travel in Z direction at the velocity v 
close to that of light, and the distribution of charge in the 
bunch do not depend on the azimuthal angle 6. 
Equations for non-zero components of plasma electrons 
momentum and electromagnetic field that describe the 
nonlinear wake-fields are (see also Refs. [2,3]): 
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As usual, Eqs. (1) and (2) were derived taking into 
account that the curl of the generalized momentum is 
zero, pH—rotP=0, or in our case 

1 He + d Pj d r-d p r l d z = 0. (7) 
Also we allow for radial variations of unperturbed 
plasma density. In Eqs. (l)-(6) y=(l+P 2 +P r 2 ) 1/2 is a 
relativistic factor, P zr =P z jY c an d N=nJn p {Q) are 
respectively dimensionless components of velocity and 
density of plasma electrons, n p (r) is the unperturbed 
density of plasma electrons, N=n p {r)ln p {Q), 
a=n b (z,r)/n p (0), n b is the concentration of bunch 
electrons, p=v/c, Vj_=3/3r+l/r. The space variables are 
normalized to X p {r=Q)l2n= Vk p (r=0), z=k p (r=0)(Z-v t), 
where X p and k p are the linear wavelength and 
wavenumber. The momenta and velocities are 
normalized respectively to mc and the velocity of light 
and the strengths of electric and magnetic fields - to the 
nonrelativistic wave-breaking field at the axis 
EJr=0)=m e (O pe (n=0)vJe (£ m [V/cm]=0.96xn/ ,2 [cm 3 ]), 
(O p =(ATtn p elm) m is the electron plasma frequency, m e 
and e are the rest mass and absolute value of electron 
charge. The field of forces acting on relativistic electrons 
in the excited field is F(-eE z , e(Hg-E), 0). 

3 THE CASE OF WIDE BUNCH 

Consider the case of wide bunch (k r b »l, where r b is the 
bunch radius) when the transverse components of an 
exciting field are small and the longitudinal components 
close to the bunch axis are approximately equal to those 
predicted by one-dimensional nonlinear theory. Here one 
can apply the perturbation method taking the ID 
nonlinear mode as the ground state. So, for wide bunches 
in uniform plasma [N p (r)=l] we shall seek the solution of 
Eqs. (l)-(6) in the vicinity of bunch axis in the form: 
P~P +pj, P~pr, E~E +l/, E~lr, H^hr, a~CC -5r, 
where P (z) and E (z) are the values of longitudinal 
momentum and the strength of electric field at the axis, 
Pzriz), /„(z), h(z), <5(z)«l. In the zero approximation in 



this values one can obtain the equations for P and E 
that describe the ID nonlinear regime that is studied 
sufficiently well [4]. In the first approximation in the 
small values we have: 

d 2 p z I d z 2 -AdpJ dz + Bp z = ft N 8, (8) 

where A=(dP/dz)2NJp, B=(Nj 'p)(pN 2 lyj-d' ' PJdz\ ft, 
N and % are respectively dimensionless velocity, density 
and the relativistic factor of plasma electrons in ID 
nonlinear regime. The remaining quantities can be 
expressed through p.. In Fig. 1 we show the focusing 
gradient fJr=(H e -E)/r=-2pJN , that is excited by d=4J 
long wide bunch with density a =0.2. 

-9t+5- 




Figure 1 : The gradient of dimensionless focusing field as 
a function of longitudinal coordinate. The bunch is d=4.7 
long, and a=0.2, 5=0.2 y={\-p 2 ) m =\Q. 

It was shown numerically based on Eq. (8) that 
although the amplitude of excited transverse components 
is small, their wavelength is nearly equal to that of one- 
dimensional nonlinear wave. The amplitude of 
oscillations periodically changes with z (see Fig. 1). The 
modulation depth grows and the modulation period 
decreases as the amplitude of nonlinear longitudinal 
oscillations increases. These effects take place also for 
other components of the field. 

4 WAKEFIELD IN UNIFORM PLASMA 

Eqs. (l)-(6) were solved numerically choosing the 
Gaussian profile of the bunch both in longitudinal and 
transverse directions: 

a = a exp[-(z- Z0 ) 2 / a 2 ]exp(- r 2 / a 2 ). (9) 

In case of small amplitudes of the excited wake wave 
(when a «l) the numerical solutions well agreed with 
linear theory predictions. Shown in Fig. 2 is the 
nonlinear 2D plasma wake wave excited in uniform 
plasma by the relativistic electron bunch with parameters 
a =0.4, (7=2, (7=5 (for example, in this case 
n M =4xl0 13 cm 3 and the characteristic longitudinal and 
transverse sizes of the bunch ojk respectively are 
1.06mm and 2.65mm when n=10 14 cm 3 ). The main 
difference here from the linear case is the change of 
shape and length of the wave with the radial coordinate 
r, as well as the change of amplitude in the longitudinal 




Figure 2: The two-dimensional nonlinear wake wave in 
uniform plasma. The parameters of the bunch are: 
a =0.4, (7=2, (7=5, y=\0. (a). The longitudinal electric 
field. 1— r=0; 2—r=A\ 3—7=6; 4— the density of bunch at 
the axis a(z, r=0). The focusing field f=H e -E r (b) and 
magnetic field strength (c). l-r=2; 2-r=4; 3-r=6. 

direction. Note also an enlargement of the range of 
focusing forces (f<0) in the nonlinear wake wave. 

It is easy to see that due to the dependence of the 
wavelength on r, the field in the radial direction grows 
more chaotic as the distance from the bunch increases. In 
fact, the oscillations of plasma for different r are 
"started" in the wake wave with nearly equal phases but 
different wavelengths. As \z\ increases, the change of 
phase in the transverse direction (for fixed z) becomes 
more and more marked. This leads to a curving of the 
phase front and to "oscillations" in the transverse 
direction. Curve 1 in Fig. 3 shows the radial behavior of 



the accelerating field E z in the nonlinear wake wave. 
Qualitatively, the radial dependence of the field differs 
from that of the linear case by the change of sign and 
"steepening" of fields along r (see also Ref. [2]). One 
can determine the longitudinal space parameter 
characterizing the wave front curving as follows: 

£ = A„/[2(1-A„/A(0))], (10) 

where A(r) is the nonlinear wavelength. At the distance 
\Az\~^ from the bunch the oscillation phase at the axis 
(r=0) is opposite to that on the periphery (r>o r ) where 
the wave is nearly linear. 
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Figure 3: The radial behavior of accelerating electric 
field strength E z . 1— E z (z=— 25, r) in the nonlinear wake 
wave excited in uniform plasma for the case given in 
Fig.2 (\AzH [see Eq. (10)]); 2-E(z=-25, r) in 
nonuniform plasma for the case given in Fig. 4. 

The magnetic field strength in the nonlinear two- 
dimensional wake wave as shown in Fig. 2(c) is different 
from zero (note that in the linear case H$=0 in the wake, 
where the wave is potential) due to the fact that the 
velocities of plasma electrons in the wave are not any 
more small in comparison with the bunch velocity. The 
magnitude of higher frequency oscillations (as compared 
to the plasma frequency) performed by the magnetic 
field along z grows in proportion to the nonlinearity. 
Such a behavior of magnetic field is a purely nonlinear 
effect. The nonlinearity of the wave implies a rise of 
higher harmonics in P z and P r . According to (7), the rise 
of magnetic field is due to these harmonics and this 
accounts for frequent oscillations seen in Fig. 2(c). On 
the other hand, according to (7), the non-zero magnetic 
field in the wake means that the motion of plasma 
electrons in the nonlinear wave is turbulent (rofP^O). 
The degree of turbulence (the measure of which is H e ) 
grows as the nonlinearity. 

From the viewpoint of acceleration and focusing of 
charged bunches in the wave, the curvature of the 
nonlinear wave front is undesired as the quality 
(emittance, monochromaticity) of the driven bunch 
worsens. Below we show that in radially-nonuniform 
plasma one can avoid the curvature of the wave front. 



5 RADIALLY-NONUNIFORM PLASMA 

Thus, in two-dimensional nonlinear regime the nonlinear 
wavelength changes with r due to nonlinear increase of 
the wavelength with wave amplitude. On the other hand, 
the linear wavelength X ~n 4 decreases with density of 
plasma. Let us assume that the nonlinear wavelength of 
the two-dimensional wake wave in the uniform plasma 
A(r) is known. Then, one can roughly compensate for 
the radial variation of the nonlinear wavelength by 
changing the unperturbed density of plasma in the radial 
direction according to the relation 

A(0) / A(r) = X p (r) I X P (0) = [n p (0)/n p (r)] 1 ' 2 . (i i) 

If we put the function A(r) to be Gaussian (according to 
numerical data for profiles (9), this is approximately the 
case at least for r«j), then one can take the transverse 
profile of the unperturbed plasma density to be also 
Gaussian: n p (r)= n p0 exp(-r /<7 2 ). It follows from Eq. (11) 
that in this case (7=r/[ln(A(0)/A(r))]" 2 . For example, the 
numerical data for A(r) in the nonlinear wave shown in 
Fig. 2 give O p ~\\. Fig. 4 illustrates the validity of this 
assertion (see also Fig. 3, curve 2). 

Numerical solutions obtained for the case of nonlinear 
wake wave excitation by a short laser pulse (Laser 
Wakefield Accelerator, see e. g. Ref. [1]) show that the 
results given in this work do not change qualitatively. 




Figure 4: Accelerating electric field of two-dimensional 
nonlinear wake wave in nonuniform plasma with a = 1 1 
for r=0, 2, 4 and 5 in the order of magnitude reduction. 
The bunch parameters are the same as those in Fig. 2. 
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